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Abstract 

We propose a new treatment of the dynamics of a periodically time-dependent Liou- 
villian by mapping it onto a time-independent problem and give a systematic expansion 
for its effective Liouvillian. In the case of a two-level system, the lowest order contribution 
is equivalent to the well-known rotating wave approximation. We extend the formalism 
to a pair of coupled two-level systems. For this pair, we find two Rabi frequencies and 
we can give parameter regimes where the leading order of the expansion is suppressed 
and higher orders become important. These results might help to investigate the inter- 
action of tunneling systems in mixed crystals by providing a tool for the analysis of echo 
experiments. 

PACS numbers: 76.60.L , 66.35.-ha 

1 Introduction 

Quantum tunneling of substitutional defect ions in alkali halide crystals leads to particular 
low temperature properties |l|, |2|. Due to their misfit in size or shape such defect ions 
are confined to a potential energy landscape with a few degenerate potential wells. At low 
temperatures thermally activated hopping is inhibited and the defect ion passes through the 
barrier by quantum tunneling. These tunneling defects are usually modeled by two-level- 
systems (TLSs), which dynamically are identical to spin degrees of freedom. 

Depending on the detailed nature of the TLS, both electromagnetic and acoustic fields might 
couple to the system. One of the interesting features of such systems is their response to 
an oscillating weak external field which is in resonance with the TLS. The key point in this 
context is the emergence of the Rabi oscillation , which is related to the amplitude of the 
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external field and which allows a systematic experimental investigation of the response of a 
TLS. 

The easiest way to derive the Rabi oscillation is to do the so called rotating wave approxima- 
tion (RWA) as originally introduced by Rabi 0]. Since then, the problem has been tackled 
by numerous methods starting with the work of Bloch and Siegert who could first give 
a correction to the RWA. Further approaches have been a solution in terms of continued 
fractions presented by Autler and Townes and especially the Floquet formalism developed 
by Shirley [^. 

In this paper, we present yet another expansion scheme for solving the problem. It is some- 
what related to the Floquet formalism, but we try to find an effective static Liouvillian or 
Hamiltonian rather than a time-evolution operator, which makes our calculations easier. Like 
the Floquet formalism, our approach can be easily generalised to systems involving several 
energy levels, where the simple geometrical interpretation of the RWA fails. Thus, we can 
define the RWA in these situations as the lowest order contribution of our expansion scheme. 

Our paper is organised as follows. In the following section, we address the time-dependent 



problem of a single spin. We first rederive the RWA, then in section 2^, we present the 
formalism for our expansion. In the third section then, we use our new method for treat- 
ing the problem of two coupled spins. The main results of this section are the emergence 
of two Rabi frequencies, whose relative magnitudes depend on the spin-spin coupling, and 
the identification of parameter regimes where higher order corrections to the RWA become 
important. 



2 The treatment of a single spin 

This section deals with the problem of one TLS that is coupled to an external force. The 
TLS is modeled by the Hamiltonian [I], || 



^sys = -k^ - , (1) 



with tunneling element k, asymmetry A, and energy e = y/k'^ + A^. The external force is 
given by 

Fext = -27? sm{^t)Y ■ (2) 

This choice of Z/cxt implies that the external field couples solely to the spatial coordinate az 
of the TLS, as realized in the case of an electric field coupling to the dipole moment of the 
TLS. 

The time evolution of the statistical operator is governed by the Bloch equation [||, ^ 

dp 



dt 



i[p,H,y,+H,,t]-l{p- p"')=:Cp , (3) 



where the right hand side of the equation defines the Liouvillian C. The second term in the 
Bloch equation describes a relaxation of the statistical operator towards its equilibrium given 
by p'^i = 6"^^=^=/ (tr e~'^^''5's j_ This choice of p^*^ is justified as long as the relaxation time 
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is much longer than the period of the driving field, 7 <C If we expand the density matrix 
in terms of the Pauli matrices and the 1-matrix, 

rl+rl + rl<l, (4) 

the Liouvillian C can be represented by a 4 x 4 matrix. Eq. (|3|) then yields a differential 
equation for the vector {l,rx,ry,rz)- 



/>=-(! + r^ax + Vyay + r^Uz 



( 1 \ 

\rz 



2.1 RWA 



Let us first rederive the RWA. To simplify the problem, we neglect the asymmetry and 
consider 

H=-e^-2r,sm{u;t)^ . (5) 
The Bloch equation in this case reads 

^=i[p,H]-j{p-p^'^) , (6) 

and p'^'' has the particularly simple form 

peg = ^(1 + nV,), = tanh(/3e/2) . (7) 

Throughout this section, we assume that both the coupling to the external field and the 
damping are small and that we are close to resonance, i.e. 

r],j,\e — u!\ ll! . (8) 



The Hamiltonian in (|5|) describes not only a tunneling system coupled to a harmonic oscil- 
lator, but equally well a spin-i particle subjected to a static magnetic field in x-direction, 
= ^/{gfJ-B), and an oscillating field in ^-direction, = 2r]/{gpB)s'm{u!t). In the latter 
interpretation, the components of the density matrix {rx,ry,rz) are proportional to the mag- 
netic moment of the spin. Thus, under rotation they transform like a pseudo- vector. In the 
following we will use the terminologies of both interpretations. 

As a preliminary, we consider the case of a static magnetic field. From the classical analogue 
one can infer that the magnetic moment will precess round the direction of the magnetic field 
with a frequency proportional to the field strength. Thus, the static Hamiltonian 

^ = ^(e.f + e,^ + e.^) , + + = 1 (9) 

describes a rotation of the vector (ri,r2,'r3) about the axis {ex,ey,ez) with an angular fre- 
quency e. The explicit formula for the statistical operator as a function of time is useful in 



3 



the following, we give it for 7 = 0, i.e. for the case without dissipation: p{t) = UQ{t)p{0) 
Q-^Ht pf^^^^iHt ^ ^j^jgpg time evolution operator UQ{t) can be represented by the matrix 



/I 







Uo{t) 



c+e^(l-c 

°y 



\ 

exey{l-c)-ezS exez{l-c) + eyS 

y 



exey{l-c) + ezS c+e^(l-c) eye^(l-c)-esS 



(10) 



\ exez{l-c)-eyS eyez{l-c)+exS c+e^(l-c) / 

using the abbreviations c := cos(et), s := sin(et). 
In the presence of dissipation, the corresponding solution reads 

Pit) = Uit)p{0) = + e-^'Uoit){piO) - p' 



(11) 



where the stationary solution p^^ is given by £ p*^* = 0. 

Now consider the Hamiltonian (|5|), which includes non-static terms. We split it into two parts 



H = 
Ho 

Hi 



Ho + Hi 
2 



-UJ ■ 



2rj sin(ci;t) 



6 :-- 



UJ 



(12) 



2 ' ' ' 2 

where assumption ^ ensures that Hi is a weak perturbation to Hq. 

The dominant part Hq of the Hamiltonian leads to a rotation of the magnetic moment about 
the X-axis. We take this into account by changing to a rotating frame of reference. In this 
rotating reference frame, we denote the statistical operator by p. It is related to p through 



Pit) 



-iHot 



p{t)e 



iHot 



p{t)e 



-iu)t(cTx/2) 



and its equation of motion reads 



7(p-p^^) =:C-p 



with the transformed Hamiltonian 



H 



j-Ta; Uy 

-0 h ?7 — 

2 ' 2 



(13) 



(14) 



(15) 



The form of the effective Hamiltonian is easily understood by observing that the linearly 
polarised field -ffext can be decomposed into two counterrotating circularly polarised fields 
of strength {Bz/2) and angular frequency u) each. In the rotating frame of reference, one of 
them becomes a static field whereas the frequency of the other one is doubled. 

In the following we want to neglect the oscillating field. This is justified by the following 
consideration. During one period of the oscillating field, p changes only little, as its typical 
time constant is l/(\/ 5"^ + rf') 3> {1/2^). Therefore, it is plausible to assume that the 
contributions of the time dependent field cancel over a period. This approximation is called 
rotating wave approximation because the original sinusoidal magnetic field is replaced by 
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a rotating field of half the amplitude and the same frequency. The remaining part of the 
Hamiltonian is static, 

^eff = -<5?^+r?5^, (16) 



describing a rotation with angular velocity (— 5, r/, 0). Thus, we can solve the equations of 
motion (14) using (^Tj) while the stationary solution reads 



P 



1 + 



n 



5"^ + rf + 7^ 



{(5^ + l'^)cFx -rjSay -Vl^^z] 



n 



tanh(/3e/2) . (17) 



If at t = 0, the system is in equilibrium, we find for the dipole moment of the TLS in the 
laboratory frame 



{p) (t) = tr ipp{t)) 



-n — 



cos 



M)(i+»- 



■It 



|2sin(17t) - 2cos(ilt)| 



+ sm{ujt) [ -4 + e~^* <{ 4 cos(ilt) + ^ sin(m) 



6j 



nn 



(18) 



using the definitions 0^ := fi'^ + 7^ := + r/^ + 7^; is called Rabi frequency. Fourier 
transformation yields peaks at uj and uj ziz Q whose heights have an overall factor of {rj/Q), 
ensuring that we only find a signal close to resonance. 

Although the method of the rotating wave approximation is intuitively very appealing, it 
lacks a rigorous foundation. In particular, there are no means to estimate the quality of the 
approximation and its limits of validity. 



2.2 An expansion scheme for the RWA 

In this section, we present a systematic expansion, which yields the rotating wave approxi- 
mation as the lowest order contribution. We consider the Hamiltonian 

F = -fc^-A^-2r?sin(a;t)^ , (19) 

where only the strength of the driving field ry and the damping constant of the Bloch equation 
7 are small compared to tu: 

ri,'y^uj. (20) 

The time evolution is dominated by the static part of the Hamiltonian. Therefore we start 
by diagonalizing this part by means of a rotation about the y-axis, 

= ~ ~ 2fr7 sin(Lijt)-^ , (21) 

where 



e = y/k"^ + A^, cj) = arctan(A;/A), u = cos(</>), v = sin((;/)) . 
After this unitary transformation, the equilibrium value of the statistical operator reads 

1 



^ 2 



1 + n°cj. 



n° = tanh(/?e/2) . (22) 
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Now, we proceed in a similar way as in the last section. The Hamiltonian is divided into two 
parts 



H = Ho + Hi 



Hq := —nu 



Hi := — (5 ^ — 2nr7 sin(a;i)^ — 2u77 sin(a;t)^, 



5:=e- 



(23) 



where the integer n has to be chosen such that \6\ < ^ thereby ensuring that Hi is small 
compared to ui. The effect of Hq can be taken into account by changing to a frame of reference 
spinning around the x-axis with the frequency \nuj\ 



Pit) 
dp 



i[p,H]-^{-p-p''i) =:Cp, 



(24) 



where L reads in the basis {l,ax,CFy,az) 



m = 

using the definitions 



/ 


V 








-7 -ur]g{t) 
wngit) -7 

—ur} f{t) —{S + 2vr) sin(a;t)) 





ut] f{t) 
6 + 2vri sin(ci;t) 

-7 



(25) 



f{t) = cos{{n — I) (jjt) — cos{{n + 1) ujt) 
g{t) = sin((n — 1) uit) — sin((n + 1) ut) . 

In the following, wc first show that there is an effective, static Liouvillc operator yielding 
a time evolution that matches the correct time evolution at integer multiples of the period 
r := I'K ju). Then, we give an expansion of the effective Liouville operator in powers of (1/cu). 

If we use the time ordering operator T, the formal solution for the evolution of p reads 



pit) = Texp 



/:{t') dt' 



piO). 



(26) 



Taking advantage of the periodicity of J0.{t), we rewrite the last equation as 



Texp 



p{t) = (rexp [ /* jC{t') dt' 
Now, we define the effective Liouvillian C^s by 

exp [Cesr] = Texp 

and end up with 



r c{t') dt' 

Jo 



N 



p{0), NT<t<{N + l)T. 



Pit) = <^ r exp 



Nt 



C{t') dt' 



r c{t') dt' 

Jo 

exp [CesNr] p{0) . 



(27) 



(28) 
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Thus, we separated the problem into two parts: first, one has to calculate the evolution 
between Nt and t accessible by simple perturbation theory. Second, one has to determine 
/^eff- We discuss the latter first. 

We use the identity 



exp[£cfTA^T] 



^ C,^NT + -Cl^{NTf + ... 

which implies that £eff is the sum of all terms on the right hand side proportional to Nr. 
Luckily, these terms form a series in (l/o;) that can be calculated explicitly order by order. 
By substituting Zi := loU, the m*^ integral yields 

2'rrN zi Zm-1 ^ 

^ J dz,J dz2 . . J dzm C{zi/io)C{z2/io) . . . Cizm/co) = ^ E Q " (27riV)^ , (30) 



where the coefficients Cj do not depend on u. Only the terms proportional to N contribute 
to £efr and we obtain 





i-Nt 


Texp 


/ C{t') dt' 




Jo 


Nt 


Nt h 


dtiC{ti)+ dh 









.(m-1) 
-eff 



1 



(31) 



at order m — 1. Our choice of 6 and the condition ( pO| ) assure that we can truncate the series 
after a few terms. 



To lowest order, we have 



Nt 



c 



(0) 
eff 



— Idt.dt,) 







This integral vanishes for n ^ 1. So, to lowest order the driving field only contributes at 
resonance (n = 1), where we obtain the results of the rotating wave approximation (see 
Section [2.11 ). A calculation of the first order corrections in the case of resonance (n = 1) 
yields 



,{0) 

-eff 



1 



^(1) 

-cfT 






^^uvP 2uvr] 



\ 



\u5 





-2uvri 




\ 



(32) 



2".^ 4"" '' ^ / 

This can be expressed in terms of an effective Hamiltonian and an effective equilibrium state 
defined by the relation 

£efrP =: i[p, H^s] - l{p - pIs) ■ (33) 

We obtain 



H. 



H. 



eff 

(0) 
eff 

(1) 
eff 



17(0) , tt{1) 
^eff +^eff 

-5- + u^- 



1 

UJ 



3 2 1 ^CTy 

—u ri uo— + Zuvr) — 

4 ' 2 2 2 '2 



(34) 
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and 



eg _ eq , (1) 
PeS — P ^ PcE 

We still have to calculate the evolution between Nt and t. This can be done using per- 
turbation theory. With Ci{t) := C{t) — Ces-, we write the time evolution operator as the 
series 



U{t) = U{t-NT)U{NT) 
t 

= e^<=ft* + y"dtie^<=«(*"*i)£i(ti)e^<=«*i 



Nt 

t ti 



+ J dtij dt2 e^-ff(*-*i)£i(ti)e^-ff(*i"*2)£i(t2)e^-ff*2 + . . . . (36) 

Nt Nt 

Thus, we see that to lowest order the time evolution is completely described by the effective 
Liouvillian, while the corrections are small due to our choice of 5 and the condition 



In the first part of this paper, we presented a formalism in which we treated a harmonically 
driven TLS. We could map the time dependent problem onto a static one and give, order by 
order, an effective Liouvillian describing it. The lowest order result could be identified with 
the RWA. 

This formalism is not restricted to a simple TLS, but can be easily generalised to any system 
driven by a periodic field. Thus we can generalise the RWA to systems that have not a 
simple geometrical interpretation like the TLS by identifying it with the lowest order effective 
Liouvillian. In the second part, we demonstrate this by investigating a coupled pair of TLSs 
in a harmonic field. 



3 A Pair of Coupled Tunneling Systems 

Depending on the detailed nature of the TLSs in mixed crystals, an electric and an elastic 
moment is connected to the defect which can thus interact with lattice vibrations, external 
fields or neighbouring defects. Consequently only at low defect concentrations one can de- 
scribe the situation by isolated TLSs. With rising concentration pairs, triples etc. of defects 
are involved until finally one faces a complicated many-body system ||2|, [l^. 

Echo experiments |^ are an experimental tool to investigate such systems [11]. In the fol- 



lowing we investigate the influence of the interaction between the tunneling defects in such 
experiments. Since we treat pairs of TLSs our results are confined to the limited range where 
the concentration is sufficiently high to give a considerable amount of pairs but still below 
the value where many-body effects set in. 

We model a pair of coupled tunneling systems (TLSs) by the Hamiltonian 

k + 5k k — 6k 

-Hpair = 7> O'a, O 1 1 (g) - JcJ^ (g) (T^ . (37) 
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The first and second term describe symmetric TLSs with k > 5k > 0, while the third term 
models the dipole-dipole-interaction of the TLSs with the coupling constant J. The external 
harmonic force, coupling to dipole moment p of the pair, is represented by 

-ffext = -r]sin{ujt){az (g> 1 + 1 03 Uz) (38) 

with the external frequency to and the coupling constant rj. This ansatz assumes that the 
two TLSs are parallel to each other and to the harmonic force. The full Hamiltonian of the 
problem reads H = //pair + -ffext and we are interested in the time-dependent dipole moment 
of the pair assuming that we start from thermal equilibrium. 

In the following, we outline the calculations, while details can be found in Appendix^ First, 
we diagonalize i/pair to find two uncoupled effective TLSs: 



-f/pair = -yo-x ® 1 - yl ® '^2;, (39) 



with the energies ei/2 = V J"^ + k'^ it + 6k'^. These effective TLSs have a more compli- 
cated coupling to the external field than the original ones, reading 

Hcxt = -V sin(u;t) {{ui + U2<7x) ®az + az(^ (us + Uiax)] (40) 

where the tt's dependend on /c, 5k and J with values in the interval [—1,1] (for further 
information see Appendix |A.1|). 



Now, we take the dominant effect of -f/pair into account by changing to a rotating frame of 
reference 

p{t) = exp{iujt{n ax 1 + ml crx)/2) p{t) exp(— iwt (n fj^, (81 1 -|- m 1 (81 Ox)l'i) (41) 
where the integers n and m have to be chosen such that 

|ei — nuj\ < -uj and |e2 — mLo\ < -lo . 

Regarding only resonance phenomena we have to distinguish three cases. The first two cases 
describe a situation where the internal frequencies ei and €2 are clearly separated. In the first 
case, the external frequency oj is in resonance with the smaller internal frequency 62 {m = 1 
and n > 1), while in the second case it is in resonance with ei (m = and n = 1). In the 
third case, lo is close to both internal frequencies (m = n = 1). This implies that both the 
energy separation 5k and the coupling J of the underlying TLSs are small compared to their 
energy k, i.e. J, 5k <C k. 

We will discuss mainly the first case with m = 1 and n > 1. The results for the second case 
are essentially the same and can be found in Appendix [A.3| . The treatment of the third case 
has also been relegated to the appendix; in spite of more complicated calculations it gives no 
principally new insight into the problem. 

Neglecting dissipation (7 = 0), focusing on the situation of exact resonance (e2 = uj) and 



following section 2.2, one obtains for the dipole moment to lowest order: 



p(0) {t) = ]y] ^Ai {sin((a; - J7,)t) - sin((u; + a)t)} (42) 
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Figure 1: The two Rabi frequencies r2i^2 foi" two different values of (5k/k) against {J/k) 



with 



1/2 



r]{ui±U2) ^1/2 = r2(l ± ri) ri/2 = tanh(/3ei/2/2) 



(43) 



Thus we find two Rabi frequencies ^^1/2 yielding four frequencies for the dipole moment: 
= £2 ± r^i and iVp = €2 i: f^2- This is one of the main results of our paper. Although 
the two TLSs can be transformed into effective uncoupled TLSs they do not give the same 
results as truly uncoupled TLSs, which have only one Rabi frequency. In Fig. |^ we have 
plotted the two Rabi frequencies (Eq. ([43|)) versus the coupling constant J. At J = 0, i.e. if 
the two spins are not coupled, the two Rabi frequencies coincide, thus reproducing the result 
of a single spin. For increasing J, l^i grows, tending quickly to 2r], whereas $^2 goes down to 
zero. In the case of negative J, l^i and 1^2 interchange their roles. 

The amplitudes of the peaks (Eq. (p3D) are proportional to the frequencies themselves, 
thereby enhancing the signal of the bigger Rabi frequency. Additionally the amplitudes also 
contain a temperature dependent factor A1/2 shown in Fig. |^. Over the whole temperature 
range, Ai is larger than A2- This is most pronounced for very low temperatures where Ai 
tends to a value of 2 whereas A2 vanishes. For higher temperatures starting at about T = ei 
they are of the same order. 

For J > 0, Ai corresponds to the higher Rabi frequency. In this case the higher Rabi fre- 
quency peak dominates the spectrum for all temperatures, reaching its maximum for small 
temperature, while the contribution of the lower Rabi frequency can be seen only at medium 
temperatures. In the strong coupling limit the lower Rabi frequency and thus also its am- 
plitude vanish, leaving us with only one Rabi frequency of 2rj. This is in accordance with 
the picture that two strongly coupled TLSs form one small effective TLS; its effective dipole 
moment is the sum of the individual dipole moments, thus yielding 0, = 2ry. 



The latter situation has been investigated by Weis et al. [11|, who measured strongly coupled 
pairs of Lithium tunneling defects in a KCl matrix. The tunneling defects, which can be 
well described as TLSs |^, IC], interact via a dipole-dipole coupling. In the special case 
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Figure 2: A1/2 = '''2(1 i n) for two different values of (e2/ei) against (T/ei) 



investigated by Weis et al. this coupling is ferromagnetic (J > 0). Due to the strong 
coupling only one Rabi peak was found. It would be desirable to investigate different pair 
constellations corresponding to weaker or even antiferromagnetic couplings in order to see 
the second Rabi frequency. 

For J < 0, Al corresponds to the lower Rabi frequency. Thus, the amplitude of the higher 
Rabi frequency is suppressed by the temperature dependent factor and we can find a param- 
eter regime, where both amplitudes are of the same order of magnitude. For strong coupling 
and low temperature, both amplitudes vanish, rendering corrections to the iiWA of order 
{r]/uj) important. These corrections lead to four new peaks close to the bigger frequency ei 
(while we are in resonance with €2!) and can be found in Appendix A. 2. Yet due to the 



overall factor of (ij/uj) we can not expect to see the new peaks in an experiment. 

In the case of resonance with the larger frequency ei, we get essentially the same results, 
but the expressions for the two Rabi frequencies and their amplitudes are different and are 
listed in Appendix |A.3| . The main difference is the existence of a parameter range in which 
we get considerable amplitudes for both Rabi frequencies. Furthermore, for |J| — > 00 both 
Rabi frequencies vanish so that again corrections to the RWA become dominant. 

Concerning the third case where both inner frequencies are close to each other, we only 
considered the situation of resonance with one of them while (ei — £2) 3> ?? (see Appendix 



A.3). The results are very similar to the previously discussed cases, but the corrections are 
of order r]/{ei — €2) rather than {rj/uj). 



4 Conclusion 



In this paper we have communicated a generalisation of the rotating wave approximation 
(RWA). First, we propose a formal framework in which for a single spin the RWA appears 
as the leading order. This formalism allows to easily calculate higher order corrections to 
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the approximation. If, on the other hand, we apply it to a general n-level-system, we can 
generalise the concept of the RWA by identifying it with the leading order contribution. 

Second we have considered a coupled spin pair in an oscillating external field. Although the 
energetic structure of such a system is well described by two isolated effective spins the time 
evolution of the system shows some 'dynamical entanglement'. These effects will disappear 
in the ferromagnetic strong coupling limit, which makes clear that our results are consistent 



with previous investigations M, O, 12 



Furthermore our calculations also show that there is a situation (concerning spin pairs) where 
terms beyond the leading order of the RWA become important. This will be the case for an 
anti-ferromagnetic coupling, where the amplitudes of the leading order are suppressed. We 
predict that whenever a situation can be modeled by a spin pair with such a coupling, rotary 
echos should show signals arising from higher order terms in the approximation scheme we 
proposed. This explicit example shows that, although the RWA seems to be correct for most 
cases, one can think of situations where the RWA does not suffice in order to separate the 
'slow part' from the complete quantum dynamics. 
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A Calculations and Corrections for the Pair 



A.l Diagonalization of the Pair Hamiltonian 

We diagonalize Hpair by the transformation 

Hpair = e-*(t'^f®'^- + ^''-®'^^) Hpa^r e^%<'v®<'-'-^''-®<'v) , (44) 

where the angles a and /3 obey 

tan(a + /?) = — t and tanfa — /?) = —— . (45) 
k ok 

This leads to the following values of u in H^xt 

ui = cos(/?) , U2 = — sin(a) , U3 = cos(q) , U4 = — sin(/3) . (46) 

A. 2 First Case: m = 1 and n > 2 

In the rotating frame of reference we calculate the effective Hamiltonian following the pro- 
cedure presented in section |2.2| . Including the dominant corrections of order |ei — u!\/u!, the 
effective Hamiltonian reads 



s s s 

Zi UJ 1 77r J. ) 



where 61 = ei — nuj and 62 = €2 — muj. Using Eq. (^6|), we can calculate the dipole moment 
up to order ij/uj. Considering exact resonance, 62 = 0, we obtain 

P^'Ht) = l{{fA^ + n^B^)sm{{u-n,)t)-{fA^-n.B^)sm{{u; + n,)t)} 
^2 + ^B,) sin((^ - n2)t) - {fA2 - IB,) sin((c. + ^2)t)} 
+ ^B2sm{ujt) (48) 
-^B^l sin((ei + r?n2)t) + sin((ei - r?n2)t) , 
-^Ssj sin((ei + r?ui)t) + sin((ei - r?ni)t) 

with 

i?i = (^^S^-^n3n4)r2(l-ri) B2 = ^{ul+ul)n = ^^u^{u:, + r2U^)n , 
while All 2 and ^1/2 as well as were already defined in Eq. (|4^). 

A. 3 Second Case: m = and n = 1 

All quantities can be obtained from those of the preceding case using the following transfor- 
mations: 

n <-> m, ri ^ r2, ui <-> and ^2 ^ tt4 (49) 

As the main results we show the Rabi frequencies plotted versus J /k in Fig.^ (corresponding 
to Fig.|l]) and the temperature dependent prefactors Ai plotted versus temperature in Fig.^ 
(corresponding to Fig.|2[). 
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Figure 3: The two Rabi frequencies ^3^4 
for two different values of {5k /k) against 
{J/k) 



Figure 4: ri(libr2) for two different values 
of (e2/ei) against (T/ei) 



A. 4 Third Case: n = m = 1 

In this case, the effective Hamiltonian reads to lowest order 

HeS = 1 - y 1 iXi cr^ + - o-y (lis + n4(T^) + - (ni + U2az) (Jy , (50) 



which is very similar to that in Eq.(^7|). Restricting ourselves to the situation where the 
energy difference of the two effective TLS's (ei — £2) is large compared to ??, we can use 
perturbation theory. If the bigger frequency is in resonance with the field, u; = ei, we obtain 
for the dipole moment: 

= i{(f^3 + 2fi?5)sin((a;-f)3)t)-(f^3-2;ii?5)sin((^ + 1^3)t)} 
- 2^^Be) sin((^ - ^4)1) - (fA^ + 2f i^e) sin((a; + ^2)1)} 

sin((e2 + VU3)t) + sin((e2 - ryua)*) 1 
-^Bgl sin((e2 + ?yM4)i) + sin((e2 - ??U4)i) 



with 



^3/4 = 'n{u3±U4) ^3/4 = ri(l lb r2) = riiiin2(l ± r2) 

B-j = {ul + ul)r2 Bf<,ioi = r2U2/i{uii2ri+U2/i) €2 = €2 + ^{ul + u^) . 



(52) 
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